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itors are Messrs. John H. Lockhart, Jr., Charles Chesnutt, and Barry W. 
Holliday. The project is under the jurisdiction of the U.S. Army Engineer 
District, Detroit. 


Work was performed under the general supervision of Ms. Joan Pope, 
Chief, Coastal Structures and Evaluation Branch (CSEB), CHL; and 
Mr. Thomas W. Richardson, Chief, Engineering Development Division, CHL. 


This report was prepared by Dr. Robert Naim and Mr. Peter Zuzek of 
W. F. Baird and Associates, Coastal Engineers, Ltd., and by Mr. Larry E. 
Parson and Dr. Andrew Morang, CHL. Field data were collected by the 
Detroit District's Grand Haven Area Office, CHL, the U.S. Geological Survey, 
Westem Michigan University, and the University of Michigan. Technical 
reviewer of the report was Mr. Edward B. Hands, CHL. 


At the time of publication of this report, Director of WES was 
Dr. Robert W. Whalin. Commander was COL Bruce K. Howard, EN. 


The contents of this report are not to be used for advertising, publication, 
or promotional purposes. Citation of trade names does not constitute an 
official endorsement or approval of the use of such commercial products. 


Chapter 1 


1 Introduction 


This report describes a study into the effectiveness of beach nourishment 
along the cohesive shore south of St. Joseph Harbor on Lake Michigan. The 
study was funded by the U.S. Army Engineer Waterways Experiment Station 
under the Monitoring Completed Navigation Projects (MCNP) Program. 


The objectives of this study were as follows: 


a. To improve understanding of the sediment transport processes for both 
fine-grain and coarse-grain sand components at this site. 


b. To improve understanding of the relationship between the movement of 
the cohesionless sediment (both fine- and coarse-grain components) and 
the irreversible downcutting of the underlying glacial till (cohesive 
sediment) at this site. 


c. To apply the improved understanding of the sediment transport and 
erosion processes in developing recommendations for beach nourish- 
ment at the St. Joseph site. 


d. To formulate general principles for beach nourishment of cohesive 
shore sites which suffer from a sediment supply deficit due to the pres- 
ence of an updrift littoral barrier. 


The study was based on a comprehensive database of the site conditions which 
were collected under the MCNP Program of the U.S. Army Corps of Engi- 
neers, Coastal Engineering Research Center. A companion report by Parson, 
Morang, and Nairn (1996) discusses geologic control on shoreline stability for 
southeast Lake Michigan. Another report by Parson and Smith (1995) 
describes an investigation of native beach characteristics for this section of the 
Lake Michigan shoreline. These supporting documents include important 
background information on the analyses and interpretations presented in this 
report, including: 


a. A description of the geologic setting. 


b. A summary of the results of the monitoring program activities. 


Introduction 


c. Laboratory experiments performed in a unidirectional flow flume to 
assess the erosion rates of undisturbed till samples which were extracted 
from the lake bed offshore of St. Joseph. 


Chapter 2 of this report provides a review of the information presented in 
the companion reports as well as an overview of the problem at St. Joseph. In 
addition, cohesive shore processes are summarized. 


A summary of the data used for the analyses completed as part of this 
investigation is presented in Chapter 3. The primary components of these data 
consist of repeated beach profiles, lake bed bathymetry, and shoreline recession 
rates. The results of additional data collection, including subsurface profiling 
with ground-penetrating radar and sediment sampling, are presented in the 
companion reports by Parson, Morang, and Nairn (1996) and by Parson and 
Smith (1995), respectively. 


Chapter 4 presents the results of a series of analyses performed to develop 
an understanding of the evolution of the shoreline and lake bed in the vicinity 
of St. Joseph and the influence of the beach nourishment program on the evo- 
lution. These analyses include 2-dimensional (2-D) and 3-dimensional (3-D) 
numerical modeling of sediment transport, profile comparisons, and bathymetry 
comparisons. 


Based on the results of the analyses described in this report, and from the 
companion report by Parson, Morang, and Nairn (1996), a descriptive model of 
the historic coastal morphodynamics in the vicinity of St. Joseph is developed 
and presented in Chapter 5. This descriptive model is used to project the 
future evolution of coastal morphology. It is in this context that the effective- 
ness of the ongoing beach nourishment program is evaluated. 


Recommendations for future nourishment efforts at St. Joseph are made on 
the basis of establishing realistic goals for the program in Chapter 6 of the 
report. A discussion of general principles for beach nourishment design on 
cohesive shores downdrift of harbor structures concludes the report. 
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Introduction 


2 Background 


Regional Coastal Processes and Geomorphology 


In general terms, this section of the southeastem Lake Michigan shore is 
characterized by eroding bluffs which consist of glacial deposits with some 
instances of relict dune formations. A detailed summary of the morphology 
and related references for this section of the Lake Michigan coast is pro- 
vided by Parson and Smith (1995). The general coastal morphology of 
Lake Michigan is described by Hands (1970). 


The lake bed also consists of glacial sediments (with isolated outcroppings 
of shale bedrock) covered with a veneer of sand and gravel of variable thick- 
ness. The sand and gravel cover represents a recent (i.e., in a geological time 
perspective) lag deposit that has been derived from the erosion of the lake bed 
and bluff in this region. Near the mouth of the St. Joseph River, the presence 
of an incised valley results in a very thick cover of sand over the underlying 
glacial sediment. However, along most of the coast, the glacial sediment is 
probably within 0 to 4 m of the lake bed surface. A discussion of the pro- 
cesses of shoreline recession on such “cohesive shores” is presented later in 
this chapter. 


The 120-year bluff recession rate, averaged for Berrien County, was about 
0.6 m/year (Hands 1976). Short-term and local rates can be much higher, 
particularly during periods of high lake levels. Downcutting of the lake bed 
between 3 and 4 m has been reported by Foster et al. (1992) for the period 
between 1945 and 1991 south of St. Joseph Harbor. The net alongshore sedi- 
ment transport direction is from north to south. The harbor jetties act as par- 
tial to full littoral transport barriers. 


Site Conditions and Beach Nourishment History 


This investigation will focus on a 12-km section of shoreline extending 
3 km north of, and 9 km south of, the harbor jetties at St. Joseph (refer to 
Figure 1). Immediately north of the harbor entrance, the fillet beach influences 
the shoreline morphology for approximately 1 km. The remaining 2 km of 
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Figure 1. Site conditions 


shoreline north of the harbor jetties are not influenced by the harbor structures. 
There is a small downdrift fillet beach immediately south of the harbor, about 
400 m in length. The 1.1-km section of shore between the fillet beach and the 
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Waterworks revetment to the south is partially protected by deteriorated groins. 
The feeder beach for the nourishment program extends from Park St. (located 
about 600 m south of the south jetty) to just south of the Waterworks revet- 
ment. Beginning with the Waterworks revetment and extending about 3.5 km 
to the south, the shore is protected by an armor stone revetment (constructed 
for the Chesapeake and Ohio Railroad over the first 1.5 km and for the high- 
way by the Michigan Department of Transportation, (MDOT), for the next 

2 km). In some places, the revetment is fronted by groins, many of which are 
in disrepair. The final 3.3 km of shore south of the end of the revetment 
consists of various forms of deteriorated wall structures and entirely 
unprotected sections. 


A Section 111 mitigation plan was implemented downdrift of St. Joseph 
Harbor in 1976 by the U.S. Army Corps of Engineers (USACE) to address the 
erosion problems that may be associated with the interception of sediment on 
the updrift side of the structures. The harbor jetties were constructed pues 
nally in 1903 and have been estimated to trap approximately 84,000 m? of 
sediment per year (USACE 1973). The mitigation consisted of placing fine 
sand from the harbor ESHA E dredging on the downdrift beaches (Johnson 
1992). More than 1,700,000 m? of sand has been placed on the beaches of 
St. Joseph. Table 1 provides the annual placement details for beach nourish- 
ment between 1970 and 1995. Parson (1992) has indicated that the fine sand 
has been a less-than-ideal material for nourishment, noting its short retention 
time and the fact that the fine sand does not fulfill the role of the coarser 
sediment which forms a large part of the natural beach closer to shore (i.e., in 
the surf and swash zones). 


Coarse material was placed on the beach in 1986, 1987, 1988, 1991, 1993 
and, most recently, in the fall of 1995 (see Table 1). This coarse sediment 
came from upland sources and was trucked to the site. The coarse grain sedi- 
ment has a ds, of about 2 mm and is well-sorted with a range of grain sizes 
from 0.1 mm to 32 mm. This material has a longer retention time and it has 
been postulated that it may protect the underlying glacial till from erosion in 
the critical nearshore zone (Parson 1992). 


The beach nourishment is placed between the ordinary high water mark 
(OHWM) (177.2 m International Great Lakes Datum (IGLD), 1985) and the 
most landward 1.2-m depth contour (174.8 m). The maximum design height 
for the placed material is 178.3 m and its eS cate width is 46 m. The typi- 
cal beach nourishment volume is about 50 m?/m over the 1-km- -long feeder 
beach or about 50,000 m? in total with fine sediment applied in the spring and 
coarse sediment in the fall. 


To classify nourishment volumes and include the results in the descriptive 
model, annual beach nourishments were grouped into three time periods: 1970 
to 1975, 1976 to 1991, and 1991 to 1995, for both fine (dredged) and coarse 
(trucked) sand (see Figure 2). Prior to the implementation of the Section 111 
plan in 1976, annual nourishment volumes averaged 23,000 m? and there was 
no trucking from inland borrow sites. From 1976 to 1991, average annual 


Chapter 2 Background 


[Table 1 
Nourishment Details for St. Joseph Harbor, Michigan 
(from U.S. Army Engineer District, Detroit) 


Dredged | Trucked Total Year 
(m3) (m3) Type (m*) Location of Beach Fill Placement 


a Ee Se ee 
cova [ Soul no. ta ee aoa] rezeonea MT pwTOR AOU J 
Ee MSE CS | 
cS AL A ll ca CA al ic 
i ae ie a CO Se | 
el Pe ca eee CN Pree | 
i nee CO ee I eee | 
El ne eo EOC eN  eree ere| 
Ce PS ACY ee es ee ee 
a ee Gee Mi Nand | 
a EE EC a ee 
[3905 | varien ana ean arose Reno PT cemioea | eae ue el 
PP Ga Un ac 
FN Ta cel eas Tl a ET 
i a Fe Cl 


1988 31 May to 29 July os 728 33, 378 Park St. - 3400 ft South 
(OHWM - 8-ft Contour) 

1988 19 Oct. to 19 Nov. 51,527 51,527 CL of Park St. - 2700 ft South 
(OHWM - 4-ft Contour) 

1989 24 May to 22 June 14,309 14,309 CL Park St. Ext 2700 ft 
(OHWM - 8-ft Contour) 

22 May to 22 June 44,515 44,515 CL Park St. Ext 2700 ft South 

(OHWM - 7-ft Contour) 

1991 3 May to 22 May 40,086 40,086 CL Park St. Ext 2700 ft South 
(OHWM 7-ft Contour) 

1991 3 Sept. to 30 Sept. 63,651 coarse 63,651 CL Park St. Ext 2800 ft South 
(OHWM - 4-ft Contour) 

1992 22 May to 9 June 25,682 fine 25,682 CL Park St. Ext 2700 ft South 
(OHWM - 7-ft Contour) 

1993 18 June to 30 July 1,756 1,756 50 ft South of CL Park St. Ext 
2700 ft S (OQHWM - 7-ft Contour) 

1993 7 Sept. to 29 Oct. 45, 821 coarse 1200 ft South of CL Park St. Ext 
1300 ft S ———— 4-ft Contour) 


Total (m° ) 1,165,374 | 588,318 
Average/year (mm?) 44,822 | 22,628 
Note: The OHWM is 1.22 m above Datum 
Denotes implementation of Section 111 Plan. 
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Chapter 2 Background 


nourishment volumes from dredging increased to 74,000 m?/year, with an 
additional 14,000 m?/year of coarse sand trucked from inland borrow sites, for 
a total of 88,000 m?/year. From 1991 to 1995, combined dredged and trucked 
volumes average only 41,000 m?/year, a substantial reduction from the annual 
volume delivered to the feeder beach between 1976 and 1991. 


Discussion of Cohesive Shores 


Sandy shores are generally distinguished by an inexhaustible local supply of 
beach sediment. In contrast, a shore is defined as cohesive when a cohesive 
sediment substratum (such as glacial till, glaciolacustrine deposits, soft rock or 
other consolidated deposits) occupies the dominant role in the change of the 
shoreline shape (i.e., through erosion). In other words, underneath any cohe- 
sionless deposit (i.e., sand and gravel) there is an erodible surface which plays 
the most important role in determining how these shorelines erode, and ulti- 
mately, how they evolve. A cohesive shore erodes and recedes because of the 
permanent removal and loss of the cohesive sediment (both from the bluff and 
the lake bed). The sand cover may come and go (depending on the season, 
water level, and storm activity), but the erosion of the cohesive layer is irre- 
versible. The characteristics of cohesive shores are discussed in more detail in 
Parson, Morang, and Nairn (1996). 


The critical point to understand is that shoreline recession, and the associ- 
ated problems of undermining of shore-based structures, could not continue 
without the ongoing downcutting of the nearshore lake bed. The long-term 
average rate at which the bluff or shoreline recedes on a cohesive shore must 
be governed by the rate at which the nearshore profile is eroded or downcut. 


Where there are downdrift erosion problems related to the interception of 
sand at an updrift barrier on a cohesive shore, downdrift mitigation efforts 
such as beach nourishment must be carefully assessed, since the sand can act 
as either protective cover or as an abrasive agent (contributing to erosion) 
depending on the quantity and type of sediment. 
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3 Existing Data Sources 


Beach Profiles 


A comprehensive profile monitoring plan was initiated with the beachfill of 
1991 under the MCNP Program. Monitoring consisted of beach profile and 
lake bed surveys taken several times each year at seven transects spaced at 
about 200 m in the immediate fill area, and additional profiles at 800-m inter- 
vals further to the south (a summary is provided in Table 2). The profile 
surveys are described by Parson, Morang, and Nairn (1996) and an associated 
sediment sampling program is presented by Parson and Smith (1995). 

Section 111 profiles are designated with the letter “R” and extend from north 
to south. Line R8 is the first profile south of the jetties, and R23 is the 
southem-most line monitored for this study (Figure 3). Four historical profile 
lines, Nos. 1, 2, 3, and 4, were also analyzed to determine multi-decade 
changes in offshore morphology. 


Lake Bed Bathymetry 


In 1995, the bathymetry of the study area was surveyed with new airbome 
technology. SHOALS (Scanning Hydrographic Operational Airbome LIDAR 
Survey) is a helicopter-mounted hydrographic surveying system which utilizes 
Light Detection and Ranging (LIDAR) to transmit and receive water surface 
and sea bottom signals. Using conventional acoustic methods, the bathymetry 
was previously surveyed in 1945/6, 1964/5, and in 1991 by the National 
Oceanic and Atmospheric Administration and the U.S. Geological Survey in a 
joint mapping project (Foster et al. 1992). 


Wave and Water Level Data 


Wave climate information has been generated by Hubertz, Driver, and 
Reinhard (1991) as part of the Coastal Engineering Research Center (CERC) 
Wave Information Studies (WIS). A detailed discussion of the wave hindcasts 
generated for this project is given in Parson, Morang, and Naim (1996). 
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Lake Michigan water levels are also discussed by Parson, Morang, and Naim 
(1996). 


Shoreline Recession 


An investigation of long-term shoreline recession rates north and south of 
the harbor jetties at St. Joseph was completed by the Land and Water Manage- 
ment Division, Michigan Department of Natural Resources (MDNR). The 
original study was completed in August 1978 to document change in shoreline 
location over a 40-year period. With the addition of a series of April 1989 
aerial photographs, the original study was updated to describe 51 years of 
shoreline change. The length of the comparison masks the influence of such 
factors as fluctuating water levels, storms, shore protection structures, and 
other natural and human disturbances. Recession data for St. Joseph are sum- 
marized in Figure 4. 


From the harbor jetties northward, the shoreline was accretional for 2.5 km, 
with an average annual accretion rate of 0.96 m/year (see Figure 4). North of 
the accretional zone, the remaining 13 km of shoreline assessed by the MDNR 
had an average annual recession rate of 0.76 m/year. South of the harbor 
jetties, only the first 0.8 km of the shoreline (corresponding to the zone of 
influence from the fillet beach) had a long-term depositional trend, while the 
remaining 13 km of shore has been eroding at varying rates. 


From the feeder beach at St. Joseph to the southem limits of the MDOT 
Revetment, the shoreline recession rates range from 0.36 m/year to 1.16 m/ 
year. There are two possible explanations for erosion along the protected 
shore south of the harbor jetties: (a) the revetment was not present for the 
entire period of the air photo comparison, and/or (b) the revetment was con- 
structed at the base of the bluff and the beach in front of the revetment has 
since eroded. When the results from the original investigation (August 1978) 
are compared to the second assessment (April 1989), in general, the annual 
rates of recession have decreased for the Railway and MDOT Revetment sec- 
tions. This suggests that the shoreline recession rate has been reduced or 
eliminated locally with the construction of the revetment. 


At Shoreham, where the shoreline is only partially protected, long-term 
recession rates are higher than to the north, ranging from 0.88 to 1.83 m/year. 
For the remaining 7 km of shoreline south of Shoreham, the average annual 
recession rate was 0.69 m/year. 
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4 Analyses of Coastal 
Processes and 
Geomorphology 


This first part of this chapter consists of three sections describing: the 
alongshore sediment transport calculations; updated cross-shore sediment trans- 
port calculations (i.e., subsequent to those presented in Parson, Morang, and 
Nairn (1996)); and the results of quasi-3-D sediment transport modeling. 


Sediment transport calculations were completed with COSMOS, which is a 
deterministic numerical model for the simulation of coastal processes. Each of 
the processes is evaluated at approximately 250 finite difference calculation 
points across the profile. The various individual predictive phases of COS- 
MOS, as well as the integrated model, have been extensively tested against 
both laboratory and field data (see Southgate and Nairn (1993), and Nairn and 
Southgate (1993)). The model is described in more detail in Parson, Morang, 
and Nairn (1996). 


The remainder of the chapter describes an investigation of the geomorphol- 
ogy of the study area through a review of nearshore profile evolution. 


Results of the Alongshore Sediment Transport 
Calculations 


Single grain size across the profile 


This section describes the results of the average annual alongshore sediment 
transport calculations that were made for each profile line using grain sizes of 
0.2 mm, 0.4 mm, and 2 mm with the original version of COSMOS-2D. Pro- 
file line locations are shown in Figure 3. 


Input for the calculation of average annual alongshore transport consists of 
a list of representative wave conditions (wave height, period, and direction) 
and durations (i.e., number of hours per year for each condition). This list was 
derived from the percent occurrence tables of height and period by direction 
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from the WIS hindcast (1956-1987). Each wave condition from the percent 
occurrence tables was run four times with the COSMOS model to represent 
four different lake levels (with the duration for each wave condition factored 
by the fraction of time associated with each of the four representative lake 
level conditions). Each input wave and water level file consisted of approxi- 
mately 1,000 conditions. 


For calculations of average annual alongshore sediment transport, the 
COSMOS-2D numerical model assumes that the profile shape remains fixed 
(i.e., profile changes due to cross-shore or alongshore sediment transport are 
not computed). The selection of input profiles is discussed in Parson, Morang, 
and Nairn (1996). 


Output from these runs consists of a description of the northerly and south- 
erly sediment transport components across each profile. Net alongshore 
transport across the profile is calculated from the two components and total 
transport for the entire profile is also calculated. Net alongshore sediment 
transport values are given for each run in Table 3. Positive sediment transport 
values represent transport to the south. All of the predicted average annual net 
sediment transport values are directed to the south. Distributions of the aver- 
age annual net alongshore transport across each profile for the three grain sizes 
are given in Figures 5, 6, and 7. While the net transport values os the 0- to 
2-mm runs fall in the range of approximately 70,000 to 80,000 m s/year 
directed toward the south, a review of the southward and northward compo- 
nents reveals that the transport is much reduced for the profiles with a revet- 
ment (i.e., R14 to R23, excluding Bee). The southward aegis transport 
component ranges from 375,328 m? at Line R8 to 170,794 m 3/year at R14. 
These differences in predicted transport rates are related directly to the profile 
shape since the same profile azimuth was assumed for each line (and since the 
same wave input was used for each line). Therefore, the low predicted values 
at Lines R14 and R23 (and to a lesser extent, R17) are a direct result of the 
deeper profiles at these locations (i.e., due to the absence of a beach at the toe 
of revetment structures). For Line R12, the peak transport occurs along the 
inner beach with a secondary peak over the first large bar. Line R14 results 
show that the peak transport occurs over the first bar offshore of the toe of the 
revetment. 


Parson, Morang, and Nairn (1996) noted that the prefill beach sediment had 
a composite ds, of about 0.3 mm and that the natural sediment (i.e., unaltered 
by beach nourishment) may be best represented by a d<, of 0.4 mm. There- 
fore, a second set of alongshore transport calculations were performed with a 
ds, of 0.4 mm. The results are summarized in Table 3 and presented in Fig- 
ure 6. For the important southward (edie transport component, te pre- 
dicted values range from 159,500 m 3/year at Line R9 to 79,900 m 3/year at 
Line Beh, This range of values corresponds more closely to the 
84,000 m 4/year which was estimated by USACE (1973) to be trapped on the 
north side of the harbor. One would expect similar values for profiles located 
north of the harbor. Sediment trapped on the north side of the harbor is 
derived entirely from the southward-directed transport component (i.e., waves 
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Table 3 
2-D COSMOS Modeling, St. Joseph Harbor, Michigan 


Average Annual Alongshore Sediment Transport 


Profile | North 


| -306,278 | 375,328 69,050 -139,824 | 154,371 14,548 -50,415 59,400 8,985 


—t 
| -291,765 | 366,579 77,814 -135,758 159,492 23,733 -44,893 2 55,421 10,528 


| -263,469 | 339,963 81,032 -122,982 | 148,291 | 25,309 -40,369 46,205 5,836 


-282,591 353,817 71,667 -131,671 | 154,095 22,424 -42,284 | 47,180 4,896 
-283,712 | 355,381 il 71,669 =| -130,454 iL 150,601 20,147 -44,679 | 48,856 4,176 
-100,086 | 170,794 -51,092 79,898 | 28,705 -16,702 24,592 7,890 
-165,290 | 243,799 -86,394 | 118,091 | 31,696 | -28,008 | 34,839 | 6,831 
| -149,609 | 255,142 105,533 -79,146 | 119,323 | 40,176 | -25,454 | 35,084 | 9,629 

-259,066 | 336,333 77,266 -124,627 | 149,683 25,506 -45,283 50,678 5,395 
| 455,01 6 | 231,740 76,724 -75,257 | 104,680 29,422 | -24,614 32,392 7,777 


from the south have little or no effect on the sediments trapped in the shadow 
of the north jetty). 


Alongshore sediment transport calculations were performed for a dsp of 
2.0 mm in the final series of these runs. The results are summarized in 
Table 3 and presented in Figure 7. Importantly, these findings indicate that 
only as little as 50 percent of the coarse sediment eroded from the feeder 
beach would make its way past Line R23 and south of the study area. The 
remaining 50 percent of the coarse sand eroded from the feeder beach is prob- 
ably deposited in the depression located offshore of the MDOT and railway 
revetments. 


Multiple grain sizes across a profile 


The COSMOS-2D model was upgraded to simulate multiple grain sizes 
across a beach profile and alongshore sediment transport calculations were 
redone for Lines R12 and R14. A dsg of 0.2 mm was assumed for the off- 
shore sediments of both profiles, with a gradual coarsening from 0.2 mm at the 
swash zone to 2.0 mm at the shoreline. 


Results for the single d,, and multiple d.q's are compared in Table 4. At 
Line R12, the COSMOS-2D model tests with a multiple grain size resulted in 
a 25-percent reduction of northerly and southerly transport from the 0.2-mm 
results of the first investigation. This is attributed to the coarsening of the 
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Figure 5. Net annual alongshore transport (ds, = 0.2 mm) 
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Figure 6. Net annual alongshore transport (ds, = 0.4 mm) 
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Figure 7. Net annual alongshore transport (deg = 2.0 mm) 
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Table 4 
Average Annual Alongshore Sediment Transport (WIS M50) 
Results from Single Grain Size and Multiple d50’s 


Average Annual Alongshore Sediment Transport 
Offshore d50 = 0.2 mm 
Beach d50 = 2.0 mm 


To the | To the Tothe | To the To the | To the To the | To the 
rte | Scam [rors _[um_[aoum_|nort [nm [ose |nort_|van_|soutt_|ncstr_|nt_| 
fara _[assaei |-2ea712| 71.66 | isos |-taoase] 20,47 [ «ease [aro [476 [essen | zener] 51 | 
[ria [17070 [-s00005|70,700 | rae | 5.002] 0,705 [2400 [16702 [7890 [1s4s04 | 97419] 7.19 | 


Note: 1. Positive transport is directed to the south. 
2. Transport calculations from average annual waves (WIS Station M50). 


sediment in the swash zone and beach for the new model runs (i.e., the coarse 
sediment has a reduced potential for transport). 


The ability of COSMOS-2D to estimate alongshore sediment transport rates 
with multiple d5,'s improves the accuracy of the predictions for St. Joseph by 
representing the natural distribution of sediment across the nearshore and beach 
zones. In general, the 0.2-mm results in Table 4 were reduced by 25 percent 
when coarse sediment was considered. For the protected sections of the 
St. Joseph shore (such as the MDOT revetment) where no beach exists, the 
reduction is less than 25 percent. 


In summary, these sediment transport calculations also indicate that perhaps 
only 50 percent of the coarse sand which is eroded from the feeder beach area 
(by storms with waves from the north) can be transported out of the study area 
south of Line 23. 


Annual variation in potential alongshore sediment transport 


To investigate variation in the wave climate, yearly estimates of wave 
energy and average direction were calculated for selected years from the WIS 
data (see Figures 8 and 9). Figure 8 shows a large annual variation in total 
wave energy ranging from a maximum in 1977 of 46,000 m/s to a minimum 
in 1986 of 17,000 m/s. The average annual wave direction presented in Fig- 
ure 9 also shows considerable variation. From the 32 years of data, seven 
individual years were selected to represent the wide range of actual combina- 
tions of wave energy and direction. 


Alongshore sediment transport rates were calculated with COSMOS-2D for 
profiles R12 (sandy shore) and R14 (revetment). Multiple grain sizes were 
considered for both profiles to represent the actual field conditions at 
St. Joseph. Results are presented in Table 5. In 1964, the average wave 
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Figure 8. Annual variability in wave energy, WIS data - Station M59 (Note: Wave energy is H2T (m?s)) 
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Figure 9. Annual variability in average wave direction, WIS data - Station M59 (Note: Shore-perpendicular azimuth at 
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Table 5 
Annual Variability in Potential Alongshore Transport (1992 Profiles) 


funnel rie 65 


To the To the To the To the 
South (m?) North (m*) Net (m*) South (m>) North (m%) Net (m?) 


Fess [zm [ouooo [arsasoareaer | -asi7 |eoraoa _|-vesseo | ssato _| 
fees [2ae _[eoec0 [easton [-aanoce | 1saoua  |aaus _|-zzrs7e _[2vasco _| 
fara [ass _[2se00 _[areare | semaa | azar __|oa7sor__|-teasc0 | e450 _| 
fe77 [2ss _[escoo [rwr.ote | seasea [aseaae [ease _[-zonsaa_[asrans _| 
fece [ace [encoo _[asaser [ease 200000 |zosoe | -eas2a [151761 _| 
ese [220 [7000 |easase [esr | sarea [vances | oreo | oose7_| 


Note: 1. Ice conditions were not considered for annual time series data. 
2. Offshore d;, = 0.2 mm and beach d,, = 2.0 mm. 
3. Wave energy (m/s) is Hs°T, (wave height squared x peak period) - provides indication of relative wave energy. 


Avg. 
Direction 


direction was 279 deg from north, which is close to the shore-perpendicular 
profile azimuth selected for the profiles at St. Joseph. Consequently, the esti- 
mated net transport for 1964 was very close to zero. Of the 7 years selected 
from the WIS data at M59, 1977 recorded the maximum wave energy 
(46,000 m/s) and the highest net southerly transport rate component of 
382,600 m?/year. 


Although the net transport rates for R12 and R14 are fairly similar, the 
southerly and northerly components at R14 are much lower than the results for 
R12, perhaps due to the deeper profile offshore of the revetment at R14. 


Historic variability in potential alongshore sediment transport 
related to profile change 


Sediment transport calculations were completed at three historical profile 
lines, Nos. 2, 3, and 4, to determine the influence of long-term profile change 
(Figure 10). The profiles were generated from the 3-D surfaces created from 
the historic bathymetry. Selection of the four profile locations was based on 
the following assumptions about the nearshore conditions and profile evolution 
prior to the comparison of the data: 


a. Line 1. Updrift cohesive profile (no influence from fillet/harbor jet- 
ties - representative of natural conditions or background erosion rate). 


b. Line 2. Updrift fillet profile (influenced by harbor jetties). 


c. Line 3. Downdrift cohesive profile (reduced sediment supply - influ- 
enced by the harbor jetties). 
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Figure 10. Location of historic profile comparisons 
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d. Line 4. Downdrift cohesive profile (not influenced by reduced sedi- 
ment supply from the north - representative of natural conditions or 
background erosion rate). 


Model runs were completed for profiles from the four available bathymetric 
surveys: 1945/46, 1964/65, 1991, and 1995. Changes in these profiles are 
discussed in a later section titled “Trends in Profile Change” (page 57), and the 
changes are illustrated in Figures 11 to 13. The results of the model runs are 
summarized in Table 6. 


Decreased depths at Line 2 in the vicinity of the depositional fillet beach 
have increased the potential for alongshore sediment transport because the zone 
of breaking waves is wider. However, this increase may be offset by the 
change in shoreline and contour orientation at this location. 


The large depression offshore of R12 and R14 and the associated steeper 
nearshore slopes at Line 3 resulted in a 24-percent reduction in southerly 
potential transport from 259,000 to 198,000 m?/year from 1945 to 1995 (see 
Table 6). At Line 3 the northerly transport component was reduced by 32 per- 
cent from 207,000 m?/year in 1945 to 142,000 m?/year in 1995. 


At Line 4, located 8.2 km south of the harbor jetties, the predicted transport 
rates were also much lower for the 1995 profile (see Table 7). From 1945 to 
1995, southerly transport decreased 31 percent and northerly transport 
decreased by 41 percent. 


In summary, long-term profile changes at St. Joseph have influenced the 
potential for northerly and southerly alongshore sediment transport. At the 
north fillet beach, the reduction in nearshore depths due to deposition has 
increased the potential for sediment transport. South of the harbor jetties, the 
deeper nearshore profiles off the revetment and the unprotected shores further 
to the south have significantly reduced the potential for northerly and southerly 
transport. 


Bypassing and channel infilling at St. Joseph Harbor 


Since the construction of the jetties in 1903, the fillet beach deposits north 
and south of the harbor have increased in size, resulting in increased potential 
for channel infilling during northerly and southerly wave attack. Profile data 
from the detailed 1995 bathymetry were used as input for the COSMOS-2D 
model to assess the existing potential for channel bypassing and/or infilling. 


Annual rates of alongshore transport were calculated beyond the end of the 
north and south jetties. Landward of the channel entrance, the harbor jetties 
due to their sheet-pile construction, were assumed to be complete barriers to 
alongshore transport. The average annual potential for channel infilling from 
the north was estimated to be 15,000 m?/year (see Figure 14a) based on the 
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Table 6 
Historic Variation in Potential Sediment Transport (Profiles from Bathymetry) 


[Wl Soeeseieeis ey coat eye cy 
To the To the Net To the To the To the To the Net 

fists [rseaso [woos [seaar lasoarr [-2orare [sates [asoaaa  [-vo7aee [sar | 
joss [2zrae [roa _|s77o1_[arsaar_|-1se7se_[secis [2saci0 _[-r7eses [57.27 | 


feet [data ro0 [esis [sasor_femnace _|resae [oes 
ses _[asose_[-vea7oe [sso [vonaso _[-aire0 [sera |rraoce [16700 [sess | 


Note: 1. Transport calculations from average annual wave conditions (WIS Station M59). 
2. Ice conditions were considered. 


1995 bathymetry and the annual average wave climate from 1956 to 1987 at 
Station M59. 


Due to lower rates of northerly transport, the south fillet beach is smaller 
than the north fillet beach. The potential rate of annual channel infilling 
during southerly wave attack is estimated to be 8,500 m?/year (see Fig- 
ure 14b). The combined annual rate of channel infilling is estimated at 
23,500 m?/year. Annual variations in wave energy could result in much higher 
channel infilling in any given year. 


The rate of natural bypassing will be significantly less than the estimated 
potential infilling rate. 


Results of Cross-Shore Modeling with Multiple 
Grain Sizes 


With the new capabilities of the COSMOS-2D model to include multiple 
grain sizes for a single profile, the cross-shore model tests for Profiles R9 and 
R14 were repeated with 0.2-mm sand offshore and 2.0-mm sand in the near- 
shore and beach. Also, an additional low water level condition was considered 
for the 24 January 1992 storm to examine the influence of low water levels on 
cohesive profile exposure and bar movement. It should be noted that the mul- 
tiple grain size version of COSMOS does not include the ability to simulate 
the mixing of grain sizes across the profile (i.e., the various grain size zones 
remain fixed in position). 


For Profile R9, located at the feeder beach, the erosion and deposition 
trends predicted with the multiple grain sizes were similar to the results for a 
single grain size as presented in Parson, Morang, and Naim (1996) (see Fig- 
ures 15 to 17). Erosion of the sand cover results in exposure of the underlying 
till in some areas. The width of the exposed till was not influenced by water 
level; however, the location of the till exposure was influenced by the different 
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Figure 14. 1995 bathymetry and transport at harbor jetties 
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Figure 15. St. Joseph, Michigan, R9 profile change for the Jan. 24, 1992 storm, low water level. 
Offshore de, = 0.2 mm, beach dey = 2.0 mm 
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Figure 16. St. Joseph, Michigan, R9 profile change for the Jan. 24, 1992 storm, actual water level. 
Offshore dz, = 0.2 mm, beach dey = 2.0 mm 
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Figure 17. St. Joseph, Michigan, RQ profile change for the Jan. 24, 1992 storm, high water level. 
Offshore de, = 0.2 mm, beach dey = 2.0 mm 
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water levels. The presence of the coarse sand (2.0 mm) on the beach protected 
this area from significant erosion, even during high water levels. However, the 
zone of transition from coarse beach sediment to fine sediment offshore was 
vulnerable to considerable erosion. 


At Line R14 the shoreline is protected by the revetment. The results of the 
profile response runs for an average water level are shown in Figure 18. For 
all water levels, there was minor erosion of the sand in the troughs between 
the nearshore bars. The trends in profile response for R14 with multiple sedi- 
ment sizes were similar to the results presented in Parson, Morang, and Nairn 
(1996) for single grain sizes across the entire profile. For all three water 
levels, the model predicted accumulation at the toe of the revetment and minor 
adjustment in the position of the large bar. 


COSMOS-3D Modeling 


Methodology 


COSMOS-3D is referred to as a quasi-3-D model because the coastal processes 
are more fully integrated across the profile than in the alongshore direction. 
The 3-D model is based on the deterministic COSMOS-2D model for the pre- 
diction of coastal processes across a nearshore profile (see Nairn (1993)). The 
profile model is extended to represent a 3-D situation through the linkage of 
11 of the individual profile lines along the St. Joseph study area. In this 
model, the profiles are treated independently in a hydrodynamic sense, but are 
linked morphodynamically by consideration of the differential rates of along- 
shore transport between adjacent profiles. Although the 3-D model grid is rec- 
tilinear, the sediment is transported alongshore in a direction coincident with 
an input “marker depth” contour. For the runs performed as part of this inves- 
tigation, two different marker depths have been utilized: (a) the most land- 
ward 2.5-m depth contour (parallel to the first large bar, which is the primary 
pathway for fine sediment transport) (b) the 0.5-m depth contour, which gives 
the alignment of the upper beach and is the major pathway for the transport of 
the 2-mm grain size. The model is only quasi-3-D and is restricted in its 
application to cases where the 3-D circulation is negligible or of secondary 
importance to the morphology change. In this respect, the St. Joseph study 
area, which features a relatively straight coast with parallel near-shore con- 
tours, is an ideal site for the application of COSMOS-3D. 


The grid for the 3-D modeling consisted of 11 of the profile lines, with 
R10a excluded due to its close proximity to R10 and R11 (Figure 19). In the 
northern beachfill area, there was 200- to 300-m spacing between profiles, 
while south of R12, the spacing was 800 m or more. The lines varied in 
length (perpendicular to the shore) between 1,000 and 1,700 m, with about 
200 calculation points describing each line. The input depth ranged from 8 to 
13 m. 
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Figure 18. St. Joseph, Michigan, R14 profile change for the Jan. 24, 1992 storm, actual water level 
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Figure 19. COSMOS 3-D setup 
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Output for COSMOS-3D consists of profile change plots for each of the 
11 profile lines, which are similar to the 2-D result format. In addition, net 
alongshore transport for the duration of the storm is plotted. The boundary 
profiles at R8 in the north and R23 in the south are assumed to remain fixed in 
the numerical calculation scheme. The results must be considered in relation 
to the large profile spacing (i.e., 200 to 800 m). The predicted changes to the 
morphology will be limited to features with lengths in the range of 200 to 800 
m or greater. Closer alongshore spacing was not possible owing to the limited 
number of survey lines. Unfortunately, the 1991 bathymetry survey was not of 
sufficient detail to pick up the shape of nearshore bars, and therefore could not 
be used to supplement the profile lines for 3-D input. The SHOALS 1995 
bathymetry was not available at the time the 3-D modeling was completed. 


Wave and water level information can be input at each of the 11 profile 
lines if information is available on variations in these parameters along the 
shore. For this investigation, the only variation that was considered consisted 
of wave sheltering effects for Lines R8 and R9 during northwest wave attack 
(i.e., in the lee of the south jetty of the harbor entrance). 


General results 


Profile change predictions from the cross-shore (2-D) modeling indicated 
that the middle and outer sections of the surf zone (which feature one or more 
bars) were relatively stable over the duration of a single storm event. A series 
of 3-D runs were performed to investigate the morphologic response of the 
sand cover under the combined influence of cross-shore and alongshore sedi- 
ment transport during storm events. These experiments were also specifically 
directed to describing the mobilization and transport of the beach nourishment 
and to assessing the exposure and downcutting of the underlying glacial till. 


Table 8 summarizes all of the 3-D results. The three initial runs consisted 
of an assessment of morphologic response under pre-fill conditions, with a 
grain size of 0.2 mm for three different storm events. These storm events are 
described in Parson, Morang, and Nairn (1996). The 2 November 1991 event 
represents one of the largest storms from the southwest (in terms of wave 
energy) over the two hindcast periods (1956 - 1987 and 1991 - 1993). The 
14 January 1992 storm featured northwest waves and was the largest storm in 
terms of wave energy over the two hindcast periods. In the first section of this 
chapter, entitled “Results of the Alongshore Sediment Transport Calculations,” 
we concluded that net transport is directed to the south owing to the predomi- 
nance of northwest storms. The 24 January 1994 event features waves which 
swung from southwest to northwest through the duration of the storm (with an 
average direction of west). This type of storm occurs frequently, and the mag- 
nitude of this particular event represents a storm that would occur once per 
year on average. The 24 January 1992 storm was used as input for all of the 
cross-shore (2-D) evaluations described in the section titled “Results of Cross- 
shore Modeling with Multiple Grain Sizes.” 
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Table 8 
Summary of 3-D Model Runs 
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Owing to the great number of output plots generated from the series of 
3-D runs (11 plots for each of the 12 storms) only the profile change results 
from run B are presented (see Figure 20). A summary of the predicted net 
alongshore transport for each of the profile lines for some of the runs is given 
in Table 9. 


Runs A to C (initial series) 


The profile change is most pronounced in the 14 January 1992 (NW) and 
2 November 1991 (SW) events. For the NW event, the alongshore transport 
values are similar to average annual alongshore transport results; alongshore 
transport is lower for the southern profiles offshore of the seawall and revet- 
ment. In other words, there is a reduction in transport moving from north to 
south, which results in deposition in the southern section, this being particu- 
larly evident at Lines R12 and R14 (see Figures 20f and 20g). In general, for 
the northwest and west storms, only 50 percent to 60 percent of the sediment 
eroded from the feeder beach area is transported beyond R23 (see Table 9 
comparing results for Lines R12 and R23). Therefore, the deep water that has 
developed through downcutting offshore of the toe of the revetment in the 
southern section of the study area acts as a partial trap to sediment moving to 
the south. 


For the southwest storms, this trend is reversed, with alongshore transport 
increasing in a northerly direction. This results in erosion between Lines R14 
and R11, which primarily affects the ephemeral beach deposit that is located 
south of the Waterworks revetment. The fact that the Waterworks revetment 
acts as a partial littoral barrier (i.e., a short groin) is not captured by the 
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a. Profile R8 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W. L. 
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b. Profile R9 3-D modeling for the 14 January, 1992 storm, dsg = 0.2 mm actual W.L. 


Figure 20. Profile change results from Run B (Sheet 1 of 6) 
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c. Profile R9a 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W.L. 
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d. Profile R10 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W.L. 


Figure 20. (Sheet 2 of 6) 
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e. Profile R11 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W.L. 
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Total Transport -108,416 m*3 


f. Profile R12 3-D modeling for the 14 January, 1992 storm, deo3 = 0.2 mm actual W.L. 


Figure 20. (Sheet 3 of 6) 
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g. Profile R14 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W.L. 
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h. Profile R17 3-D modeling for the 14 January, 1992 storm, de, = 0.2 mm actual W.L. 


Figure 20. (Sheet 4 of 6) 
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Total Transport -73,101 m*3 


i. Profile R20 3-D modeling for the 14 January, 1992 storm, dg = 0.2 mm actual W.L. 
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Total Transport -58,739 m*3 


j. Profile R22 3-D modeling for the 14 January, 1992 storm, deg = 0.2 mm actual W.L. 


Figure 20. (Sheet 5 of 6) 
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k. Profile R23 3-D modeling for the 14 January, 1993 storm, de, = 0.2 mm actual W.L. 


Figure 20. (Sheet 6 of 6) 


3-D model results. Sand eroded from the beach south of the Waterworks 
revetment (and from the south end of the feeder beach) is deposited at the 
fillet beach. 


The predicted change for the 24 January 1992 event is not that much differ- 
ent from the changes predicted under the 2-D modeling, which is a result of 
the relatively low net alongshore transport values associated with this storm, 
which swings from southwest to northwest. However, erosion in the southern 
section of the feeder beach is predicted with increasing southerly transport in 
this area (from Line R10 to R11). This erosion is balanced by some minor 
deposition in the vicinity of Lines R12 to R14. 


The increased “volatility” of the sand cover that occurred in the NW and 
SW storms results in both more and less exposure of glacial till compared to 
the 2-D cross-shore results of the section titled “Results of Cross-shore 
Modeling with Multiple Grain Sizes” (i.e., where only the inner surf zone fea- 
tured significant changes to profile shape). With the 3-D results, larger areas 
of till were exposed where a section of the profile was subject to erosion due 
to increasing alongshore transport, whereas depositional conditions occurred in 
other areas, burying till that was previously exposed in troughs between bars. 
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Table 9 


3-D COSMOS Modeling, St. Joseph Harbor, Michigan! 


ae Water Level Actual Water Level High Water Level Actual Water Level 
dsp = 0.2 mm deo = 2.0 mm deg = 0.2 mm dsp = 2.0 mm 
Marker Depth = -2.5 m Marker Depth = -2.5 m | Marker Depth = -2.5 m | Marker Depth = -0.5 m 


Pre-fill J | Pre-fillK | Post-fill L 

2-Nov-91 | 14-Jan-92 | 24-Jan-92 | 2-Nov-91 14-Jan-92 | 2-Nov-91 14-Jan-92 | 14-Jan-92 14-Jan-92 
Pa freass | vere [aaa] 2ase [sess [ezoor | essa [rezee [aces 
re fraser | asear [seat este [ease fasos | rieeo [7.200 [076 
rea Jouve [ress [rae | sae0 [ease [anros[rosow [ean aria 
fo fone | aera [ve0e | se00 [arse —[eaore | aos0r aaio [reo 
fa Jensoe | ener fosss [cao [sane |eaars | nao [verze [070 
ao Jessie | raror [azo | oes [asarfrasoe | nase [ra00e [2000 
ze fanase | saree [roa asee [iste [esr | soese [re ber 
zs [soaor | enase | veze0 [ear [raore Jerars | rare | 005 [vas] 
ee 


! Negative is southward for 3-D analyses. 
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Runs D to F (post beachfill series) 


The three initial runs were repeated with Lines R9 to R11 augmented with 
50 m?/m of fine sand (0.2 mm) beach nourishment. This represents a low to 
average level of beach nourishment with a total volume of about 50,000 m? 
extending from about 2.4 m above datum to 1.2 m below datum. The results 
generally featured very little change in the predicted alongshore transport rates 
compared to the initial prefill series. It may be recalled that the orientation of 
the nearshore contours in the numerical model are established based on the 
2.5-m contour (i.e., the “marker depth”). Therefore, since the beach nourish- 
ment only extends down to a depth of 1.2 m, in the numerical model experi- 
ments, the beach nourishment does not have an influence on the rate of 
alongshore transport related to the orientation of the contours. Nevertheless, 
for the NW waves of the 14 January storm, slightly greater deposition is pre- 
dicted at Lines R12 and R14 than in the initial series. Cross-shore transport 
processes result in rapid profile adjustment with erosion of the upper beach 
and deposition offshore for each of the storm events. For storms with 
NW waves, the erosion of the feeder beach does not appear to result in the 
rapid movement of a defined pulse or slug of sediment. 
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Runs G and H (2 mm, coarse sediment) 


The 3-D version of COSMOS used here requires a single representative 
grain size. The first two series of runs (A-F) were completed with a grain size 
of 0.2 mm. The SW and NW storm events were then repeated with 2-mm 
sediment. This coarser grain size is representative of the coarse beachfill 
derived from upland sources. However, the coarse fraction is generally found 
only close to shore (see Parson and Smith (1995)) and therefore, the results of 
these runs are probably only valid for the inner surf zone and upper beach 
areas. 


The reduction in predicted alongshore transport (compared to the fine grain 
size) was not as dramatic in these runs as it was for the average annual along- 
shore transport results (see Table 9). In general, the profiles are more stable, 
as expected. Therefore, areas of exposed till remain exposed, while sections of 
buried till are not uncovered. 


Runs | and J (high water level) 


In this series of runs, the 3-D change for the SW and NW storms is pre- 
dicted under high lake level conditions (the grain size of 0.2 mm represented 
pre-fill conditions). For the SW storm event, the predicted alongshore trans- 
port is lower for the southern section from Line R14 to R23 than that predicted 
with the average lake level conditions (see Table 9). The alongshore transport 
is reduced by 10 to 50 percent owing to the greater water depths offshore of 
the revetment. In contrast, the alongshore transport rates for the northern 
section of the study area shoreline (including the feeder beach and the fillet 
beach) are increased by about 40 percent as larger waves can reach steep sec- 
tions of the upper beach under the high lake level conditions. For the SW 
storm, these changes to the alongshore transport result in greater erosion south 
of the Waterworks revetment (i.e., Line R12) and at the north end of the 
feeder beach (R11), which is balanced by greater deposition along the north 
feeder beach and fillet beach. 


These trends are reversed for the NW storm results, with increased 
southward-directed transport in the northern section and slightly decreased 
southerly transport in the southern section of the study area. As a result, pre- 
dicted deposition is slightly greater at Lines R12 and R20. 


Runs K and L (pre- and post-fill with 0.5-m marker depth) 


As noted earlier, placement of the beachfill did not change the contours 
which the incident waves encounter for 3-D Runs D to F. This was a result of 
the fact that the “marker depth” which delineates the main pathway for along- 
shore transport was specified as the 2.5-m depth contour, which is located well 
offshore of the toe of the beachfill. Therefore, Runs K and L were performed 
with a marker depth of only 0.5 m, which better represented the contour 
changes created by the beachfill. Evaluations of pre- and post-beachfill 
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conditions were performed with a grain size of 2 mm. The coarse grain size 
was selected because the 0.5-m marker depth represents conditions where most 
of the sediment is moving relatively close to shore. 


Review of the modeling results (summarized in Table 9) shows that slight 
increases in alongshore sediment transport are predicted in the beachfill area 
for the post-fill Run L (Lines R9a to R11) compared to the pre-fill Run K. 
However, most of the profile change in the beachfill area for Run L is related 
to profile adjustment. The gradients in potential alongshore transport are very 
low and do not result in rapid redistribution of the feeder beach sediment to 
the south. 


Results of Runs H and K, which differ only in the assigned marker depth 
(2.5 m versus 0.5 m, respectively) are very similar, with only minor deviations 
in the beachfill area (see Table 9). 


Assessment of predicted glacial till downcutting 


For each of the 3-D model runs, downcutting of the exposed areas of 
glacial till was determined based on the magnitude of shear stress and the rate 
of wave energy dissipation at the location of the exposures, over the duration 
that the till was exposed. The approach used to calculate downcutting is 
described in Parson, Morang, and Nairn (1996). 


Downcutting results for the 2 November 1991 (SW) storm are compared in 
Figures 21a to 21k, showing the difference in predicted downcutting between 
the 0.2-mm and 2-mm results (i.e., Run A versus G). In general, the predicted 
downcutting (or vertical erosion) of the exposed till areas is predicted to be in 
the range of 0 to 0.15 m. The occurrences and magnitude of downcutting are 
greater for the 0.2-mm sediment for the northern section of the shoreline (1.e., 
where a beach exists with till underneath). For the southern profiles offshore 
of the revetment, the magnitude and occurrence of downcutting are similar. 
This finding relates to the fact that changes to the sand cover were much more 
limited for the lake bed offshore of the revetment for both fine- and coarse- 
grain sediment. Modeling was also performed for the 14 January 1992 storm, 
and the results were similar to the 2 November 1991 results. 


The downcutting predicted for these storm events corresponds to storm 
conditions that might be expected once every | to 5 years. The bathymetry 
comparisons have identified annual lake bed lowering in the range of 0.06 to 
0.10 m/year for the revetment shoreline between 1945/6 and 1964/5. Although 
lowering rates decreased between 1964/5 and 1991, lowering in the range of 
0.09 to 0.13 m/year occurred between 1991 and 1995. The numerical model 
results of this investigation imply that significant downcutting is still ongoing, 
not only in the vicinity of Profiles R9 to R12, but also for the deeper profiles 
offshore of the revetment between Lines R12 and R23. 
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a. Profile R23 cohesive profile erosion during 2 November, 1991 storm, dg, = 0.2 mm & 2.0 mm 
actual W.L. 
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b. Profile R22 cohesive profile erosion during 2 November, 1991 storm, deg = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. Downcutting results for the 2 November 1991 (SW) storm, Runs A and G (Sheet 1 of 6) 
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c. Profile R20 cohesive profile erosion during 2 November, 1991 storm, de, = 0.2 mm & 2.0 mm 
actual W.L. 
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d. Profile R17 cohesive profile erosion during 2 November, 1991 storm, de, = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. (Sheet 2 of 6) 
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e. Profile R14 cohesive profile erosion during 2 November, 1991 storm, deg = 0.2 mm & 2.0 mm 
actual W.L. . 
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f. Profile R12 cohesive profile erosion during 2 November, 1991 storm, d.. = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. (Sheet 3 of 6) 
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g. Profile R11 cohesive profile erosion during 2 November, 1991 storm, de, = 0.2 mm & 2.0 mm 
actual W.L. 
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h. Profile R10 cohesive profile erosion during 2 November, 1991 storm, deg = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. (Sheet 4 of 6) 
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i. Profile R9a cohesive profile erosion during 2 November, 1991 storm, dg, = 0.2 mm & 2.0 mm 
actual W.L. 
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j. Profile R9 cohesive profile erosion during 2 November, 1991 storm, deg = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. (Sheet 5 of 6) 


55 


Chapter 4 Analyses of Coastal Processes and Geomorphology 


depth below datum (m) erosion of cohesive profile (m) 
[0] 
— Initial Profile 


— Cohesive Profile 

~~ Final (dS50=0.2mm) 
—-Coh. Erosion (0.2mm) 
~ 7 Final (dS50=2.0mm) 
—- Coh Erosion (2.0mm) 


-0.2 
600 800 1,000 1,200 1,400 


distance (m) 


k. Profile R8 cohesive profile erosion during 2 November, 1991 storm, de, = 0.2 mm & 2.0 mm 
actual W.L. 


Figure 21. (Sheet 6 of 6) 


The findings also indicate that the 2-mm grain size sediment is much more 
effective than the 0.2-mm grain size sediment at mitigating downcutting for 
those sections of shore where a beach deposit protects an underlying till layer. 
For profiles which feature deep water at the toe of a shore protection structure 
(and no beach), the coarse grain size sediment is no more effective than the 
fine grain size sediment in protecting the underlying glacial till from 
downcutting. 


Summary of the 3-D results 


A major limitation of the 3-D modeling was the limited number of profiles 
available to describe the bathymetry along the study area shore. The 3-D 
modelling was completed prior to the availability of detailed bathymetry from 
the 1995 SHOALS survey. The implication of this limitation was that the 
model results could only be interpreted in a general manner; detailed changes 
to the bathymetry were either not predicted or were not entirely reliable. Not- 
withstanding this limitation, several conclusions can be made based on the 
results of the 3-D experiments. 


The 3-D results confirmed that the deep water located offshore of the 
southern revetment-protected section of shore creates an impediment to along- 
shore transport. For northwest wave events, this results in some deposition in 
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this area, including deposition at Line R12 which relates to the extension of 
the ephemeral beach south of the Waterworks revetment. 


Southwesterly storms resulted in redistribution of sediment from the beach 
immediately south of the Waterworks revetment and from the south end of the 
feeder beach northwards to the fillet beach area. 


The beachfill was predicted to respond very rapidly in the cross-shore 
direction, with the upper beach sediment being eroded and transported offshore 
as the profile readjusted to equilibrium form. However, the influence of along- 
shore transport had less immediate effects on the redistribution of sediment 
outside the feeder beach area. The sediment was predicted to move along the 
shore at a relatively slow rate, with deposition only perceptible at Lines R12 
and R14 during the NW storm event. This result is in part related to the issue 
of actual versus potential alongshore transport rates for Lines R8 and R9. If 
the actual transport rates are Jess than the potential values predicted for 
Lines R8 and R9 during a northwest storm event, the volume of sediment 
transported to the feeder beach from the north will be much lower than the 
model predictions. A reduction in the rate of sediment transport from the 
north will accelerate erosion of the feeder beach and the associated alongshore 
transport of sediment may occur more rapidly than predicted during a north- 
west storm event. 


The influence of alongshore transport on the movement of bars can result in 
significant changes to the exposure of glacial till over the duration of a single 
storm. Presumably, it is these changes which contribute to the ongoing down- 
cutting of the underlying glacial till in the vicinity of thick bar deposits. The 
volatility of the sand cover is diminished along the southern section of shore, 
which features deeper water offshore of the toe of the revetment. Predicted 
downcutting rates are much lower in this area and only occur for isolated 
sections of lake bed. The 3-D runs with a 2-mm grain size support the 2-D 
findings: because of the stability of the sand cover, existing exposures of till 
remained exposed and buried sections remained protected. The influence of 
fluctuating lake levels on the exposure of the underlying glacial till has also 
been shown to be an important factor (see section titled ““Bathymetry Compari- 
sons and Sediment Budget Calculations”). 


Trends in Profile Change 


In the past 50 years, several factors have influenced the volume of sand 
above the cohesive profile at St. Joseph, including: obstructions to alongshore 
sediment transport (harbor jetties), construction of shore protection structures, 
the Section 111 beach nourishment program, and annual variability in along- 
shore sediment transport. The quantity and stability of the sand cover above 
the glacial till has an important impact on the magnitude and location of cohe- 
sive downcutting. Long-term profile comparisons were made from the four 
snapshots of the lake bed bathymetry (1945/6, 1964/5, 1991, and 1995) and are 
discussed below. A review of the profile data collected from 1991 to 1995 
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was completed by CERC for the shore south of the harbor jetties to assess 
short-term trends in the changing volume of sand and gravel above the 
cohesive profile. 


Long-term profile change 


Four locations were selected for the long-term comparison of beach profiles 
(see Figure 10). 


At almost 3 km north of the St. Joseph harbor jetties, Line 1 is located 
outside the zone of influence of the harbor structures. From 1945 to 1965, 
there was severe lowering of the lake bed (see Figure 22). This lowering was 
probably the result of erosion of the underlying till. In contrast, the 1965 to 
1995 comparison showed little or no lowering but does feature bar migration 
related to water level fluctuations and wave action. 


Line 2 is located 1,350 m north of the harbor jetties, in a transition zone 
between shoreline influenced by the fillet beach and unaffected shoreline. Due 
to the low density of soundings in the 1945/6 survey, the bar features cannot 
be discerned. However, a depositional trend is evident from 1945 to 1995, 
especially in the nearshore zone (see Figure 11). The water level at the time 
of the 1995 survey was approximately 1 m above the historic lows recorded in 
1964/65. The 1995 nearshore sandbar is located approximately 100 m inshore 
of the 1964/5 nearshore bar (see Figure 11). 


Line 3 is located in the middle of the large offshore depression in the lake 
bed that developed between 1945 and 1995. Severe erosion (vertical displace- 
ment) of the sand and cohesive substrate occurred between the 10-m depth 
contour and the shoreline (see Figure 12). The 5-m depth contour moved 
inshore by 450 m in 50 years, for an average annual contour recession rate of 
9 m/year. With the construction of the revetment several decades ago, the 
shoreline position has been fixed at Line 3; however, the nearshore profile has 
continued to erode and the nearshore slopes have become progressively steeper 
from 1945 to 1995. 


Line 4 is located 8.2 km south of the jetties in a zone which, up until 
recently, may not have been significantly influenced by the harbor jetties. 
From 1945 to 1965, profile lowering occurred from the shoreline out to the 
-7-m depth contour (see Figure 13). There appears to be some recovery of the 
nearshore profile from 1965 to 1991, which corresponds to the period when 
the Section 111 beach nourishment program was introduced updrift of Line 4. 
A second possible explanation for the apparent gain in nearshore sand levels is 
the onshore migration of bar features due to water level rise between 1964/65 
and 1991 (refer to Parson et al. (1996)). From 1991 to 1995, the profile com- 
parison revealed significant lowering of the nearshore profile over a 
200-m-wide zone. 


Chapter 4 Analyses of Coastal Processes and Geomorphology 


0 


0 


8 


(ul) s0u¥jSIG 


S66} 0 9b/SPE} ‘UOSHedWOD ajyo1d oUO\SIY | eUIT “Zz ounBI4 


=] 
oO 
s 
> 
ee] 
oO 
ry 
= 
o 
= 
= 
=} 
So 


Sr6l 


59 


Chapter 4 Analyses of Coastal Processes and Geomorphology 


60 


Short-term profile change 


Changes in profile volumes were calculated by CHL for 14 stations moni- 
tored over a 4-year period from 1991 to 1995 (Figure 4). Three distinct zones 
for the volume calculations were selected: beach/nearshore bar, offshore bar, 
and offshore. The three zones identified distinct characteristics of the offshore 
zone south of St. Joseph. However, given the size of the study area and the 
diversity of the nearshore conditions, the locations of these zones could not be 
standardized by distance (shore-perpendicular) or depth. Consequently, their 
location on each profile was based on individual morphology. In all cases, the 
earliest profile was used as a baseline to compare changes in the volumes of 
sediment on the profile from 1991 to 1995. Positive volumetric changes indi- 
cate the amount of sand on the profile has increased since the base year of the 
comparison. Negative volume changes occur when there is a reduction in sand 
cover and/or irreversible lowering of the cohesive profile. For the wide off- 
shore zone in particular, it is noted that small errors in the profile surveys 
could result in large errors in estimated profile volume changes. 


An overlay of the long (beach and offshore) profiles for R8 is presented in 
Figure 23. Line R8 is located in a transition zone between the south fillet 
beach and feeder beach. The bathymetry comparison showed that the lake bed 
north of R8 was stable or accretional. This area is in the lee of the harbor 
jetties and any sand transported into this zone is effectively trapped (1.e., 
because of sheltering from northerly wave attack). 


The results of the profile volume calculations and the timing and volume of 
beach nourishment for Profile R8 are presented in Figure 23. The beach/ 
nearshore bar showed deposition from 1991 to 1995; however, the offshore bar 
continued to erode despite the beach nourishment. 


Profile R10 (Figure 24) is located in the feeder beach zone. The volumetric 
results suggest that the beach nourishment has been successful in maintaining 
the profile volumes in all three zones: beach/nearshore bar, offshore bar, and 
offshore (see Figure 24). Over 100,000 m? of sand was placed on the feeder 
beach in 1991, which initially added approximately 300 m?/m of sand to the 
R10 profile. In 1993, the volume of sand on the profile increased to 
600 m?/ m above the volumes recorded on August 14, 1991 (the base profile). 


Profile R11 (Figure 25) is located at the southern end of the feeder beach 
and is in a transitional zone between depositional and erosional profiles. All 
three profile zones, especially the offshore bar, experienced an erosional trend 
from 1991 to 1995 (see Figure 25). One possible explanation is the influence 
of the nearshore beach slope, which increased from 1:85 at R8 (which was 
depositional in the beach/nearshore bar zone) to 1:30 at profile R11. The 
numerical modeling indicated significant quantities of sediment could be trans- 
ported in a cross-shore direction during a severe storm, especially when the 
nearshore slopes are steep. 
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During the 3-year period (1991 to 1994) of the comparison at ie (Fig- 
ure 26), the volume of the beach/nearshore pay increased by 160 m 3/m, while 
the offshore bar volume increased by 500 m 37m. The bathymetry comparisons 
revealed that R14 is located in the zone of most severe lake bed erosion, 
exceeding 4 m from 1945 to 1995. However, during the period of the profile 
volume comparison, R14 may have been a depositional sink for sediment from 
the feeder eee (see Figure 26). If the annual rate of accumulation (approxi- 
mately 200 m 3/m) at R14 had occurred over a 200-m Segue of the shore, the 
depression in the lake bed near R14 could trap 40,000 m? of sand per year 
moving in an alongshore direction, which is approximately equivalent to the 
annual volume of beach nourishment from 1991 to 1995. 


Profile R22, located 5 km south of the harbor jetties, also appears to be 
directly influenced by the MOUS AeA on the feeder beach (Figure 27). In the 
fall of 1992, after 130,000 m? of beach nourishment was placed on the feeder 
beach in 1991 and early 1992, the SSeS ite eS pat and offshore bar gained 
significant quantities of sand, 375 m 37m and 220 m?/m, respectively (see Fig- 
ure 27). From the fall of 1993 to the fall of 1994, the trend reversed and the 
beach/nearshore bar and offshore bar eroded below the base volume of 
August 30, 1991. This erosion trend may be explained by a break in the nour- 
ishment program from the spring of 1992 to the fall of 1993, decreasing the 
rate of sediment available for alongshore transport to the beaches south of the 
harbor jetties. 


Exposure of the Cohesive Substrate 


Exposure and downcutting of the cohesive profile underneath the sand or 
gravel lag at St. Joseph are the fundamental processes that determine at what 
rate the shoreline retreats over time. Several factors can lead to the exposure, 
or increase the potential for exposure, of the cohesive profile: (a) water level 
fluctuations and associated bar migration in response to wave action, (b) reduc- 
tion in the overlying sand/gravel cover, (c) increase in nearshore beach slopes, 
and (d) changes in sediment grain size. The latter three factors have been 
investigated in the previous sections of Chapter 4, and in general, they have 
failed to fully explain the relatively even distribution of nearshore downcutting 
evident from the bathymetry and profile comparisons. This even distribution 
requires that the underlying till is exposed at all locations at some time. Vari- 
ations in bar position with changing water levels appear to provide the missing 
explanation. 


Monthly and yearly fluctuations in mean water levels for Lake Michigan 
are described in Figure 28. During the period of the investigation, 1945 to 
1995, there was extreme variability in lake levels, with a low yearly mean of 
0.3 m below chart datum (IGLD ‘85) recorded in 1964 and the high of 1.3 m 
above chart datum in 1986. On Great Lakes shores, rising lake levels together 
with wave action move the bar formations onshore and conversely, during 
falling lake levels, bars move offshore. Continuous migration of the bar and 
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trough features in the nearshore zone in response to fluctuations of water levels 
distributes the exposure and downcutting of the glacial till across the entire 
profile. 


In a technical paper prepared for the U.S. Army Corps of Engineers, Hands 
(1979) compared profiles compared for a 55-km stretch of shoreline on eastern 
Lake Michigan. From 1967 to 1976, the profiles clearly showed the shoreward 
migration of bar formations with rising lake levels (see Figure 29). Therefore, 
changes in the position of bars, and the troughs between the bars where the till 
is often exposed, result from changes in water levels. The range of water level 
variation on Lake Michigan, explains how the downcutting can be distributed 
across the entire shoreface. 


Bathymetry Comparisons and Sediment Budget 
Calculations 


The authors recognize that comparisons of bathymetry that was mapped for 
navigation purposes can sometimes produce misleading results due to the rela- 
tive inaccuracy of these surveys. However, the extent of lake bed change at 
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Figure 29. Nearshore bar migration, eastern Lake Michigan (Hands 1979) 
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this site is such that this is not an issue. Lake bed lowering is on the order of 
meters, which is far greater than any possible relative datum or measurement 
errors in the surveys. 


From the extent of the hydrographic surveys, eight panels (or sediment 
compartments) were created north and south of the harbor jetties. Wherever 
possible, the panel boundaries were selected to delineate changes in the lake 
bed evolution characteristics (i.e., south fillet beach, Panel 3 and feeder beach, 
Panel 4). Data from historic and recent bathymetries were used to create 3-D 
surfaces representative of the lake bed conditions at the time of the surveys. 
From the comparison of the 3-D lake bed surfaces, net volume changes were 
calculated for each individual panel (i.e., a positive value is obtained when the 
net change for a panel was deposition for the period of the comparison). To 
provide a relative basis of comparison for the calculated change in panel vol- 
umes between the individual periods (i.e., 1945/56 to 1964/65 and 1964/65 to 
1991) and panels, the total change in panel volume is divided by the surface 
area of the individual panel. This provides an averaged depth change for the 
entire panel, representing an annual rate of erosion or deposition. 


1945/46 to 1964/65 


The 1945 to 1965 bathymetry comparison provides a description of the lake 
bed evolution downdrift of the harbor jetties before the implementation of the 
Section 111 Plan for beach nourishment at St. Joseph (Figure 30). 


North of the jetties, the average annual lake bed lowering for Panel 1 was 
2.7 cm/year. This rate of lake bed lowering is very similar to Panel 8 (2.6 cm/ 
year), at the southern limit of the bathymetry comparison (see Table 7). This 
may indicate that Panels 1 and 8 were representative of the background rate of 
erosion at St. Joseph from 1945 to 1965. Panel 2 corresponds to the northern 
fillet deposit, adjacent to the northern jetty at the mouth of the St. Joseph 
River. The average increase in lake bed elevation was almost 2 cm/year, with 
1.0 - 3.0 m of deposition recorded near the northern end of the jetty over the 
20-year period. 


The southern fillet deposit, Panel 3, experienced minor volume increases, 
amounting to an average of 0.5 cm/year for the panel, with the majority of the 
deposition located adjacent to the shore and southern harbor jetty (see 
Figure 30). 


From Panels 4 to 7, severe lowering of the lake bed occurred between 1945 
and 1965. From profiles R8 to R17, a 3-km-long depression in the lake bed 
developed, with lowering in excess of 4 m recorded over the 20-year period 
(see Figure 30). Average lake bed lowering for the four panels (4 to 7) ranged 
from 5.7 - 10.0 cm/year. 


The total average annual sediment volume lost from Panels 3 to 8 was 
258,000 m?/year between 1945 and 1964/65. Assuming the long-term average 
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Figure 30. St. Joseph, Lake Michigan contours of lake bed elevation change, 1945/46 to 1964/65 
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annual net transport rate of 80,000 m?/year to the south is blocked by the 
jetties, and must therefore be eroded from the lake bed and shoreline south of 
the jetties, this leaves an additional 178,000 m?/year of lake bed erosion that 
must be related to offshore losses. It is likely that the majority of these losses 
are related to offshore dispersal of silt and clay associated with the erosion of 
the cohesive sediment. From these simplified sediment budget calculations, 
the NOES annual erosion of the cohesive profile may have been as high as 
20 m 3/m/year for this 9-km-long section of shoreline south of the harbor 
(Panels 3 to 8). 


1964/65 to 1991 


The second period of lake bed evolution corresponds with the introduction 
of beach nourishment at St. Joseph. From the lake bed surface change plot 
given in Figure 31 and the volume estimates in Table 7, the 1964 to 1991 
comparison clearly represents a period of reduced erosion rates combined with 
nearshore deposition, both north and south of the harbor jetties. There are 
several possible explanations for the nearshore deposition: (a) migration of 
nearshore bars during the high water levels in 1991; (b) deposition associated 
with the alongshore transport of beach nourishment; and (c) error in bathyme- 
try or datum conversion. 


Since the 1991 bathymetry did not extend north of the jetties, the 1964/65 
bathymetry was compared to the 1995 survey for Panels 1 and 2. Although 
Panel 1 continued to erode, the average annual lowering was only 12 percent 
of the 1945 to 1965 rate (see Table 7). The depositional trend also continued 
for Panel 2; however, only at 27 percent of the 1945 to 1965 rate. The long 
narrow depositional feature in Panel 2 (see Figure 31) is associated with a 
change in bar location (profile change) in response to the difference between 
the 1964/5 low water levels and the average water levels in 1995. 


At the southwestern corner of Panel 2, which corresponds with the end of 
the north jetty, the 1964/65 to 1991 bathymetry comparison in Figure 31 
revealed localized deposition in the range of 1 to 2 m for the 27-year period. 
A similar trend was also evident in the 1945/6 to 1964/5 comparison, although 
the zone of high deposition was located closer to the shore (see Figure 30). A 
decrease in the offshore depths at the end of the jetty is the direct result of the 
growth of the fillet beach deposit. This process may also have contributed to 
the development of a sediment pathway for channel infilling during north- 
westerly wave attack. 


For Panel 3 south of the jetties, the annual erosion rate averaged 0.2 cm/ 
year. Panels 4 through 7, which experienced the most extreme erosion 
between 1945 and 1965, continued their lowering trend, but at dramatically 
reduced rates. The annual lowering rates for Panels 4 to 7 ranged from 0.2 cm 
to 1.1 cm. The total lake bed lowering for the isolated case of the large 
depression in the lake bed (R12 to R20) was in the range of 1 - 2 m for this 
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period. The zone of the greatest lowering rates also shifted southward from 
R9 - R17 to R12 - R20 (refer to Figures 30 and 31). 


A comparison of Figures 30 and 31 clearly shows that the highest rates of 
deposition in the nearshore between 1965 and 1991 correspond to the locations 
of significant lowering between 1945 and 1965. This would suggest that sand 
moving in a cross-shore and alongshore direction is partially trapped in the 
deeper nearshore zones, until the depressions are filled. Consequently, these 
depressions may be sinks for sand transported from the feeder beach. 


The net volume change related to the changes in the lake bed surface 
(within the panel boundaries) south of the harbor jetties was -25,000 m?/year 
between 1964/65 and 1991. Disregarding the negligible sediment input from 
shoreline recession and channel SESS Uy input to the sediment budget from 
the beach nrouals nga averaged 88,000 m 3/year from 1965 to 1991. Assum- 
ing that 80,000 m $/year of the nourishment sediment is transported in an 
alongshore direction (due 1} the southerly directed net transport gradient), and 
that approximately 8,000 m 3/year may have been lost as annual deposition in 
the navigation channel, there would be no net gain or loss resulting from 
alongshore transport processes. Therefore, the annual net OREM 1 in the sedi- 
ment budget from the bathymetry comparisons of -25,000 m 3/year must be 
largely related to offshore losses. As noted above, offshore losses are probably 
the result of the erosion of the cohesive sediment and the offshore dispersal of 
silts and clays. Consequently, for the panels south of the harbor jetties, the 
volume Qe irreversible erosion of the cohesive substrate may have been as high 
as 2.8 m 3/m/year. 


1991 to 1995 


The 1991 to 1995 bathymetric comparison is limited to Panels 3 to 8 south 
of the harbor jetties, as seen in Figure 28, due to the limited surveying done in 
1991. During these 4 years the volume changes south of the jetties changed 
dramatically. Volumetric losses and lake bed lowering were greater than the 
previous peak during the initial interval (1945 to 1965) (see Table 7). 


The most dramatic erosion rates between 1991 and 1995 occurred in a 
200-m-wide band along the shoreline, with 1-4 m of lake bed lowering. Depo- 
sitional areas, seen in Figure 32, are further offshore and do not compensate 
for the nearshore lowering. There are several possible explanations for the 
reduced rate of offshore deposition: (a) sand eroded from the nearshore is 
widely dispersed offshore; (b) sand eroded from the nearshore is transported in 
an alongshore direction; (c) a significant percentage of the eroded nearshore 
volumes is glacial till and provides very little sand to the local sediment 
budget. 


The volume loss related to lake bed lowering was approximately 
367,000 m /year for this period. Disregarding the negligible inputs from 
shoreline recession (outside of calculated panel volumes) and harbor 
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Figure 32. St. Joseph, Lake Michigan contours of lake bed elevation change, 1991 to 1995 
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bypassing, approximately 80,000 m? annually must be eroded from the lake 
bed to supply the potential for net southerly alongshore transport and an 
additional 8,000 m? is lost annually to deposition in the navigation channel. 
Considering es the annual nourishment velinmies between 1991 and 1995 
were 41,000 m?/ year, approximately 47,000 m 3/year of lake bed erosion 
would be required to supply the additional losses i alongshore transport and 
channel infilling. Consequently, of the 367,000- m? loss related to lake bed 
lowering, up to 320,000 mm may have been attributed to pe irreversible 
lowering of the cohesive profile or approximately 35.5 m 3/m/year. 


1945/46 to 1995 


Figure 33 compares 1945/6 and 1995 bathymetry, and represents 50 years 
of lake bed evolution at St. Joseph. With the exception of the northern and 
southern fillets (Panels 2 and 3), the entire lake bed has experienced dramatic 
lowering. A large depression in the lake bed has been created by 2 to 5 m of 
vertical erosion in the nearshore zone between St. Joseph and Shoreham. 
Given the approximate size of the depression, over 3,000,000 m> of sediment 
has been eroded from the lake bed in the last 50 years. 
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Figure 33. St. Joseph, Lake Michigan contours of lake bed elevation change, 1945/46 to 1995 
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Interpretation of Results - 
A Descriptive Model of 
Coastal Morphodynamics 


This chapter provides a summary of the study findings in the form of a 
description of the historic, present, and possible future coastal processes and 
morphologic evolution for the study area shoreline and Iake bed. Two periods 
are considered for the historic changes spanning the 1945/6, 1964/5, and 1991 
bathymetry surveys. The present conditions are represented by the changes 
between 1991 and 1995. The study area shore may be subdivided into seven 
sectors as follows (see Figure 34): 


a. 


A section of coast north of the harbor which appears to be uninfluenced 
by the presence of the harbor (i.e., this corresponds to Panel 1 of the 
lake bed surface comparison analysis). 


The updrift fillet beach located immediately north of the harbor jetties 
(i.e., Panel 2). 


The downdrift fillet beach extending about 400 m south of the harbor 
(i.e., Panel 3). 


The feeder beach area extending from Line R8 to the Waterworks 
revetment (i.e., Panel 4). 


A section with uninterrupted shore protection in the form of revetment 
and seawall from the Waterworks revetment to Line R22 (i.e., Panels 5 
and 6). 


The section of unprotected or partially protected shore extending from 
Line R22 to south of Shoreham; (i.e., Panel 7). 


A section at the southerly limit of the study area which, historically, 
does not appear to have been influenced by the harbor jetties (i.e., 
Panel 8). é - 
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The descriptive model of coastal morphodynamics for the historic, present, 
and future conditions at St. Joseph is illustrated in Figures 35 - 37. For each 
period, the descriptive model summarizes the findings of the lake bed surface 
comparisons, dredging and nourishment volumes, recession rates, and repre- 
sentative estimates of alongshore sediment transport. Any variation in the 
transport rates between the different periods only reflects changes to the repre- 
sentative nearshore profiles used as model input with the average annual wave 
climate from 1956 to 1987. Qualitative information from the analysis of lake 
levels, profile comparisons, and COSMOS 2D/3D modeling presented in Chap- 
ter 4 is also incorporated in the descriptive model. 


Historic Conditions 


1945/6 to 1964/5 


Descriptive model results for 1945 to 1965 are summarized in Figure 35. 
Current understanding of coastal processes for cohesive environments would 
suggest that the shoreline of Sector A has been eroding since the glaciers 
receded several thousand years ago. Analysis of air photo information dating 
back to 1939 indicates that the recent long-term recession rate is about 0.8 m/ 
year (see Section of Chapter 2 entitled “Shoreline Recession”). The average 
lowering rate for Panel 1 was found to be 2.7 cm/year for the period from 
1945 to 1965. This rate of erosion compares well to the situation in Sector G 
during this period (i.e. at the south end of the project area), which featured a 
recession rate of about 0.9 m/year and an averaged lowering rate of 2.6 cm/ 
year. Based on this finding, and on the fact that these recession rates are 
similar to those found in areas further to the north and south of the harbor 
(i.e., well beyond any zone of harbor influence), the authors suggest that these 


two sectors are representative of the “background” erosion conditions related to - 


cohesive shore processes and are not strongly influenced by the presence of the 
harbor jetties, at least for historic periods. This is an important finding 
because the Section 111 program is only intended to mitigate erosion related to 
the presence of the structure and not the background erosion. 


The fillet beaches immediately north and south of the St. Joseph River 
mouth (Sectors B and C) have been stable or accreting at least since the con- 
struction of the jetties in 1903. Numerical modeling results indicate that sig- 
nificant quantities of sediment may be deposited in these areas during storms. 
The bypassing analysis showed that the combination of the long jetties and the 
deep navigation channel acts as a total littoral barrier, trapping all sediment 
reaching this area from either the north or the south. 


Somewhere in Sector D (i.e., the feeder beach), the shore changes from 
sandy to cohesive as the bank of the incised river valley is encountered. 
MDNR calculations of long-term recession rates indicate that the entire reach 
of Sector D has been eroding, with recession rates between 0.35 and 1.15 m/ 
year over the last 50 years (with the larger rates occurring immediately north 
of the Waterworks revetment). 


Chapter 5 Interpretation of Results - A Descriptive Model of Coastal Morphodynamics 


79 


LEGEND 


Modeled Longshore Transport 
Potential in m3/yr (COSMOS 2D/3D) 


Deposition Outside of 


Deposition in Panel #2 


Bathymetry Depression 


From Bathymetry 
Comparisons 


cam = Accretion + 
Erosion 


Dredging .: Se , 


Accretion During Southerly ———-=— 
Wave Attack o” 


Growth of Fillet 
Outside of Panel 


-1.16 m/yr 


LAKE q : No Nourishment 
MICHIGAN 


ST. JOSEPH 


Large Depression In Lake Bed 
Developed Between 1945 — 1965 


-0.46 to 
Erosion of Overlying Sand -0.98 m/yr Annual Recession 
Vaneer and Erosion Ly Rates (Mich. DNR) 
of Glacial Till 


Rapid Lake Bed Lowering 
(2 to 3m/20yr) 


Panel #8 Lowering Similar 
to Panel #1 (representative of R23 
background erosion rate) SHOREHAM 


vi © 0.79 t 
) mm -0. (-) 
¥ as rey -1.83 m/yr MAP SCALE 
cy Fee 1:50,000 
ae 1000m 500m 0 
aif ES 


. Growth of Fillet Deposits (Panels 2 & 3) 6. 1964 Water Level was all time low 

. High Variability in Long Term Water Levels 7. No Bypassing, only Channel Infilling 
. Lowering in Panels #1 & #8 is similar 8. Transport Volumes are m3/yr 

. Panels 4-7 severe lowering (downdrift deff.) 

. Cohesive Lowering Estimated at 20 m3/m/yr south of harbor jetties 


Figure 35. 1945 to 1965 descriptive model 
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Figure 36. 1965 to 1991 descriptive model 
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Figure 37. 1991 to 1995 descriptive model 
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Sector E consists of a 4-km-long stretch of uninterrupted shoreline protec- 
tion. There is a transition zone between the feeder beach and the deeper pro- 
files offshore of the revetment protecting the railway further to the south. The 
profile for Line R12 shows that significant lake bed lowering in front of the 
revetment at this area has not yet occurred (see Figure 30). However, further 
offshore, 4 to 5 m of lake bed lowering occurred from 1945 to 1995. Both the 
numerical modeling results and observations from aerial photos show that the 
beach located immediately south of the Waterworks revetment is subject to 
large fluctuations. The revetment itself probably acts as a groin structure 
impounding at least some sediment to protect the underlying glacial till in the 
nearshore zone most of the time at Line R12. 


Sectors F and G consist of the section of coast extending from Line R22 to 
south of Shoreham. Here, the shore is only partly protected or entirely unpro- 
tected. This section features long-term recession rates of 1 to 2 m/year 
between 1945/6 and 1964/5. 


1964/65 to 1991 


The 1964/65 to 1991 period is characterized by much lower rates of deposi- 
tion or erosion in the nearshore zone (see Figure 29) when compared to the 
earlier 1945 to 1964/5 period (see Figure 30). The possibility that the 1991 
bathymetry featured an error in vertical control or datum conversion was inves- 
tigated and dismissed as a possible explanation for the discrepancy between the 
rates of change between the two periods. An extensive review of all original 
data and datum conversions applied to the hydrographic surveys did not iden- 
tify any errors. The observation of low erosion rates in Sectors A and G, 
which were previously identified as representative of the background erosion 
condition, coupled with a low deposition rate in the Sector B fillet and erosion 
in the Sector C fillet located south of the harbor, suggests that the driving 
force of erosion and deposition (i.e., wave-driven sediment transport) may have 
been less effective than during the previous period (see Figure 36). Unfortu- 
nately, the available wave climate information only extends back to 1956, and 
it is not possible to substantiate this hypothesis. 


A more certain explanation for reduced lake bed lowering rates in Sectors 
D, E, and F is the influence of the Section 111 beach nourishment program, 
which was initiated in 1976 (with some nourishment placed as early as 1970). 
In these sectors, there was a tenfold decrease in the lake bed lowering rates. 
In Sector G, representative of background conditions, the lake bed erosion rate 
was lower by a factor of only 2.5. The trend for this period suggested that the 
Section 111 Program was successful in mitigating the lake bed lowering rates 
for Sectors D to F. While it may be argued that a beneficial effect was also 
experienced in Sector G, it is more likely that the reduced erosion rate in this 
sector can be explained by generally lower driving forces during this period as 
mentioned in the previous paragraph. 
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It should be noted that the tenfold reduction in lake bed lowering rates may 
not be sufficient with respect to mitigation of the harbor influence on erosion 
further downdrift. If the feeder beach sand simply ends up slowly filling the 
large lake bed depression that has developed in Sectors D to F, the shore fur- 
ther downdrift will continue to be denied the historic levels of sediment 


supply. 


Existing and Future Conditions 


Existing conditions (1991 to 1995) 


Comparison of the lake bed surfaces from the 1991 and 1995 hydrographic 
surveys reveals a rapid acceleration in lake bed lowering. In Sectors E to F, 
the lake bed lowering rates are 30 to 50 percent higher than the 1945/6 to 
1964/5 comparison period and an order of magnitude greater than the 1965 to 
1991 period (see Table 7). Of greatest concem is the observation that the rate 
of lake bed lowering in Sector G is of a similar magnitude to that of Sectors E 
to F (see Figure 37). In other words, it would appear that Sector G is now 
being influenced by the harbor structure and may no longer be regarded as 
representative of background erosion. A review of the contour plots of lake 
bed change (see Figure 32) also indicates that the focus of lake bed lowering 
(i.e., that led to the development of the depression offshore of Sectors E to F) 
has shifted to the south. 


One significant difference between this most recent period and the previous 
comparison period was the annual average volume of beach nourishment. 
Annual placement volumes have been reduced by approximately 50 percent to 
40,000 m? over the last 5 years (see Figure 2). The reduced level of beach 
feeding may at least partly explain the accelerated erosion rates. 


Projections of future conditions 


The fillet beach south of the harbor is currently stable or slightly accreting. 
During southwest storms, this sector receives sediment from erosion in the 
feeder beach area. It would appear that the fillet has reached its maximum 
extent and that any additional sand transported northwards eventually makes its 
way into the navigation channel where it is deposited, and later dredged. 


The feeder beach shoreline is maintained at a stable average position with 
the annual beach nourishment. Without the nourishment, the numerical model 
investigations have shown that shoreline recession would recommence, with 
the transport of sand to the south and the uncovering and downcutting of 
underlying glacial till where it exists. Comparison of the 1991 to 1995 lake 
bed surfaces in Figure 32 revealed that this sector experienced erosion under 
the recently reduced nourishment levels. A summary of the changes to the 
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profile lines in this sector (i.e., Lines R8 to R11) over the period from 1991 to 
1993 is given in Table 10. This table indicates that each beach nourishment is 
followed by rapid profile adjustment (PA) or moderate to high erosion (ME to 
HE). 


Table 10 at) 
St. Joseph Harbor, Lake Michigan, Summary of Profile Data (Beach Fills, Profile 
Change, and Wave Energy) 


Profile influenced by beach fill i LD low depostion 

first profile survey i Sl MD = medium deposition 
no change between profiles i i HD high deposition 
Profile adjustment (net change is zero) 


Wave W.L. Above | Volume 
Energy Datum (m) 


(m2 5)* | (GLD‘8S) | (m*3/m) 


Note: Wave energy (m’s) is H? * T (wave height squared x wave period) provides indication of relative wave energy. 
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The transition area of Sector E from the Waterworks revetment to south of 
Line R12 is distinguished by an ephemeral beach feature (i.e., a beach subject 
to significant fluctuation in size). Table 10 indicates that in two of the three 
years of the monitoring program, this sector experienced deposition towards 
early fall (as a result of the beachfill moving south) and erosion in late fall as 
the deposit was eroded by subsequent storms. Results of the model tests 
indicate that this sector is subject to the highest alongshore transport rates for 
the study area shoreline. Numerical model results also indicate that this area is 
subject to ongoing downcutting, particularly offshore of the beach deposit. As 
noted in the section entitled “1945/6 to 1964/5,” the Waterworks revetment 
probably helps to impound sediment and maintain the beach immediately south 
of the revetment. Numerical model tests also revealed that the coarse-grained 
beachfill derived from upland sources is much more effective at protecting the 
glacial till under the beach in this sector. 


Although downcutting of the nearshore profile in Sectors E and F may 
eventually diminish owing to the deeper water that has developed offshore of 
the shore protection, the numerical model results suggest it is still ongoing, as 
did the 1991 to 1995 lake bed comparison (see Figure 32). Model results also 
indicate that there may be ongoing deposition of sand in this sector, since only 
about 50 percent of the coarse sediment transported into this sector from the 
north is predicted to be transported southwards beyond Line R23. It was seen 
that during the 1965 to 1991 period with higher annual beach nourishment 
volumes, the rates of nearshore profile lowering in this sector were signifi- 
cantly reduced (see Table 7). Table 10 shows that this sector typically 
receives sediment sometime in mid to late fall. Based on the predicted reduc- 
tion in potential transport rates of coarse sediment between Line R12 and Line 
R23, we estimate that about half of the 600,000 m? of coarse fill that has been 
placed since the beginning of the Section 111 nourishment program has been 
deposited in this sector. With this assumption, and assuming the deposition 
occurs over a 500-m-wide band of the shore extending out to the 6-m contour, 
the average gain in thickness of sand cover would be 0.067 m since 1976. 
Based on the findings of lake bed surface comparisons and the results of the 
numerical model tests, this annual deposition rate of 0.0035 m/year derived 
from the beach nourishment is at least balanced, and probably outpaced by the 
ongoing downcutting of the underlying glacial till. This was certainly the case 
during the 1991 to 1995 period, with lower annual beach nourishment 
volumes. 


In order to raise the profiles to the historic lake bed levels (i.e., to allow 
unimpeded sediment transport to the south), and assuming about half of the 
traditional coarse beach nourishment volume (i.e., about 20,000 m’/year since 
1986) is deposited in this sector and that downcutting can be arrested in the 
near future, almost 8 million m’ of sediment would be required over the next 
400 years at the current rate of nourishment. The numerical model tests indi- 
cated that the 2-mm grain size sediment was no more effective than the 
0.2-mm sediment in protecting the underlying till from exposure and down- 
cutting in this sector. 
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The most southerly sector, Sector G, and the unprotected shoreline further 
to the south received about half of the historic net alongshore sediment supply 
rate of coarse sediment. This is because the deep water offshore of Sector D 
acts as a sink for about 50 percent of the coarse beachfill sediment. Therefore, 
it is likely that the shoreline of Sector G, and particularly the shore to the 
south of this sector, are suffering due to a depletion of the historic sand cover 
with the associated increased exposure of the underlying till and increased 
downcutting and shoreline recession rates. The loss of the coarse fraction 
results in greater erosion close to shore (i.e., where slopes are steeper and only 
the coarse-grain-size fractions remain relatively stable under most conditions). 
The most recent lake bed comparison (1991 to 1995) revealed that the lower- 
ing had in fact increased dramatically in Sector G compared to earlier periods. 


Comments on the Effectiveness of the Beach 
Nourishment Program 


The fillet beach of Sector C would probably remain stable without beach 
nourishment from the Section 111 program. At present, perhaps as much as 
50 percent of the sand placed in the feeder beach area (particularly for the 
dredged finer sediment) ends up back in the navigation channel from where it 
was originally removed (and will be removed again). 


There must be a more cost-effective approach to maintaining the position of 
the shoreline in Sector D than beach renourishment. An alternative approach 
may also be more environmentally acceptable and less disruptive to the local 
community (i.e., not requiring the annual trucking operation for the placement 
of coarse sand and gravel). 


The primary local beneficiary of the ongoing nourishment is the transitional 
part of Sector E. Here, too, there may be more cost-effective means of pro- 
tecting this section of shoreline. The coarse sediment is much more effective 
than the fine in protecting the till underneath the beach in this sector. The 
coarse sediment fulfills a role (which would have been present historically) in 
protecting the underlying till from downcutting that the fine sediment cannot 
(i.e., over the steeper nearshore slopes). 


Sector E has been a sink, possibly for up to 50 percent of the coarse sedi- 
ment placed in the feeder beach area. However, the effectiveness of this 
sediment (whether the coarse grain or fine grain type) in counteracting the 
ongoing downcutting (either presently or in the future) is questionable. There 
may be more cost-effective means of protecting the toe of the existing struc- 
tures. It is unlikely that the placement of the 8 million m* of beach nourish- 
ment required to completely. fill the depression that has developed over time is 
justifiable. 
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During the period from 1986 to 1995, Sector G and the area to the south 
received perhaps 50 percent of the coarse sediment eroded from the feeder 
beach. Therefore, this sector and the shoreline to the south experience a 
deficit compared to the historic sediment supply. This situation, combined 
with the depleted supply during the years prior to 1976, must have resulted in 
decreased sediment cover in this area and may have caused an increase in 
downcutting and shoreline recession. Comparison of the 1991 and 1995 lake 
bed bathymetries indicates the problem of accelerated offshore lowering and 
the related shoreline recession has extended south of Sector G. 


It would be much more effective to place the entire annual allotment of 
beach nourishment (or at least the trucked coarse sediment) south of Lines R22 
or R23 where it would be 100 percent effective in supplying the downdrift 
shores. The erosion problems in the study area could be addressed with site- 
specific solutions. With this action, the implementation of further shoreline 
structures to the south of Line 22, to counteract the increased erosion, may be 
avoided. 


Recommendations for Future Monitoring 


The following monitoring activities should be continued to assess the effective- 
ness of modifications to the beach nourishment program. 


a. Aerial photos should be continued to monitor the level of shoreline pro- 
tection in and south of Sectors F and G. 


b. Aerial photos should be regularly analyzed to monitor recession rates in 
and south of Sections F and G to update the MDNR data. 


c. Lines R12 to R23 and new lines further to the south should be moni- 
tored regularly to improve understanding of the lake bed changes in 
these areas. 


d. A complete survey of the lake bed, both north and south of the harbor 
jetties, should be completed 5 to 10 years after the 1995 SHOALS 
survey, or after significant modification to the beach nourishment 
program. 
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6 Beach Nourishment Design 
Guidelines 


Based on the findings of this investigation and the knowledge of cohesive 
shores that has developed since the early 1980’s, some general design 
guidelines are presented for the specific circumstances of St. Joseph, and for 
some general categories of cohesive shore situations. 


Recommendations for St. Joseph 


Lowering of the lake bed offshore of the MDOT and C&O revetment (i.e., 
Sector E in Figure 34) is a result of both the interruption of alongshore 
transport (particularly prior to the initiation of the Section 111 program) and 
the stabilization of the shoreline position related to the construction of the 
revetment. 


The present beach nourishment program does not appear to provide any 
significant benefit to the stability of the revetment along the Sector E shoreline 
or to the lake bed offshore of the revetment. This is despite the fact that 
perhaps 50 percent of the beachfill sediment is deposited permanently on the 
lake bed in this sector, and volume losses dropped to less than one fifth their 
former 20-year average during the 30 years after nourishment was initiated. 


Beach nourishment is definitely effective at maintaining a stable shoreline 
position in Sector D. The coarse grain sediment is an essential component 
which protects the till under the upper beach from downcutting during storms. 
Fine-grain nourishment on its own (i.e., from dredging alone) is, however, 
insufficient to protect the underlying till from exposure and downcutting. 


Placement of unrestricted beachfill (i.e., without any substantial retaining 
structures such as headlands) is probably not a cost-effective means of main- 
taining an average stable shoreline position. A solution to retaining a 
permanent beach at this location should be sought through the use of rock 
headlands or breakwaters. It may be argued that this is not the intention of the 
Section 111 program; however, it must be recognized that this has been the 
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result of and would continue to be the result of an unmodified nourishment 
program. 


The greatest flaw in the current nourishment program is that the area where 
a supply of sediment is most urgently required is only receiving 50 percent or 
less of the historic supply rate of coarse sediment. This seems to have 
accelerated recession rates for the shoreline south of the study area (i.e., 
Sectors G and southward 1991 - 1995). These erosion pressures result in 
construction of more shoreline protection by property owners. In the long 
term, these actions only further aggravate the problem by further reducing the 
supply rate (by eliminating the input of sediment from shoreline erosion and 
by impeding alongshore transport as deep water develops offshore of the struc- 
tures). 


The authors recommend that beach nourishment be placed downdrift of 
Line R22 so that 100 percent of the fill reaches the area where it is required 
(i.e., versus the current situation where perhaps 50 percent or less of the coarse 
beach nourishment is deposited in Sector E without any apparent benefits). 
The nourishment should consist of both fine (dredged) and coarse grain 
components. By moving the feeder beach to the south, the sedimentation rate 
experienced in the navigation channel should be significantly reduced. As a 
result, maintenance dredging costs may be reduced if less frequent channel 
dredging is needed. 


General Recommendation for Beach Nourishment 
on Cohesive Shores Downdrift of Harbor 
Structures 


It must be recognized that cohesive shores have very different erosion 
characteristics from sandy shores and this has a significant impact on the 
downdrift nourishment requirements. In addition, there are varying degrees of 
cohesive shores (related to the extent and role of the overlying sand cover), 
which also have an important influence on the nourishment requirements. 


Furthermore, effective downdrift nourishment requirements must be 
determined in light of changes to the lake bed that may have occurred as a 
result of the presence of the harbor structures prior to the initiation of a 
nourishment program. This is not necessarily the case for sandy shores down- 
drift of harbor structures. 


Beach nourishment guidelines for the two extremes of cohesive shore 
conditions (with respect to extent of historic, predevelopment sand cover) are 
discussed here. A final special condition is also considered. 


In some cases, sections of cohesive shore on the Great Lakes (and else- 
where) will feature only a “limited” sand cover. As a possible defining 
variable, the sand cover between the 4-m depth contour and the bluff would 
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have a volume of less than 100 m?/m in these cases. Under these conditions, 
the underlying glacial till is either only thinly covered (i.e., with beach and bar 
thickness of less than 1 m) or entirely exposed. In other words, the till is 
frequently exposed over the entire profile to conditions of active downcutting. 
In these situations, it is not clear that the impoundment of sand in an updrift 
fillet beach, and the deprivation of this sand from the downdrift beaches and 
lake bed will have any measurable impact on the rate of lake bed downcutting 
and the associated rate of shoreline recession. This hypothesis was 
successfully applied in the Port Burwell (north central shore of Lake Erie) 
litigation case where the Govemment of Canada successfully defended against 
a $30-million claim which held that the harbor structures at Port Burwell had 
caused accelerated recession for 40 km of downdrift cohesive shore (see 
Philpott (1986)). 


The opposite extreme consists of a situation where the glacial till under- 
neath the sand cover is rarely, if ever, exposed in the natural condition (prior 
to the construction of harbor jetties). This situation has been documented for 
the Illinois shoreline north of Chicago by Shabica and Pranschke (1994). In 
this case, the interception and impoundment of alongshore sediment by large 
shore-perpendicular structures has resulted in a reduction of sand cover from 
over 500 m?/m to less than 200 m?/m in places. In this case, the reduced sand 
cover resulting from the impoundment at the shore-perpendicular structures 
results in accelerated shoreline recession along the downdrift shore. Beach 
nourishment is required in these cases, not only to reinstate the historic 
sediment supply rate, but also to replenish the sand cover to its historic level. 
The latter requirement may be achieved through augmenting the sand cover 
volume to its natural level (this may not be practical or realistic owing to the 
large volumes required). Otherwise, the requirement may be relaxed if the 
effectiveness of the protective characteristics of the overlying sand cover can 
be augmented. The protectiveness of the sand cover could be improved 
through the provision of sediment which is coarser than the natural or native 
sediment. Specific grain size requirements should be determined based on the 
profile shape, properties of the underlying till, wave exposure, and sediment 
transport characteristics (both alongshore and crossshore). 


A special condition of cohesive shore which may be relatively common 
relates to cases where the natural profile shape is convex instead of concave 
(see Stewart and Pope (1993)). Gray and Wilkinson (1979) document the 
existence of this type of cohesive shore at locations on the east shoreline of 
Lake Michigan north of St. Joseph. This condition is a result of the presence 
of a more erosion-resistant surface in the nearshore. The protected nearshore 
shelf may consist of some form of bedrock or glacial till that is armored by a 
boulder and cobble lag deposit. Shoreline (or bluff) recession on this type of 
cohesive shore is particularly sensitive to changes in lake level. While 
downdrift nourishment requirements for this type of cohesive shore may be 
less in volume (i.e., less than what might be determined based on potential 
transport rates), the timing and grain size characteristic requirements should be 
carefully considered. 
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In summary, the nourishment requirements for cohesive shores downdrift of 
harbor structures (or other impediments to alongshore transport) are more 
complicated than the requirements for similar situations on sandy shores. The 
requirements must be established on a site-specific basis. They may vary from 
cases where no beach nourishment is required to others where the natural 
supply must be completely replaced and/or augmented with coarse grain 
sediment. 
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